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ABSTRACT 



Context. Ground-based imaging and imaging spectropolarimetric data are often subjected to post-facto reconstruction techniques to 
improve the spatial resolution. 

Aims. We test the effects of reconstruction techniques on two-dimensional data to determine the best approach to improve our data. 
Methods. We obtained an 1-hour time-series of spectropolarimetric data in the Fe i line at 630.25 nm with the Gottingen Fabry-Perot 
Interferometer (FPI) that are accompanied by imaging data in the blue continuum at 431.3 nm and Can H at 396.85 nm. We apply both 
speckle and (multi-object) multi-frame blind deconvolution ((MO)MFBD) techniques. We use the "Gottingen" speckle and speckle 
deconvolution codes and the MOMFBD code in the implementation of Van Noort et al. (2005). We compare the resulting spatial 
resolution and investigate the impact of the image reconstruction on spectral characteristics of the Gottingen FPI data. 
Results. The speckle reconstruction and MFBD perform similar for our imaging data with nearly identical intensity contrasts. MFBD 
provides a better and more homogeneous spatial resolution at the shortest wavelength when applied to a series of image bursts. The 
MOMFBD and speckle deconvolution of the intensity spectra lead to similar results, but our choice of settings for the MOMFBD 
yields an intensity contrast smaller by about 2 % at a comparable spatial resolution. None of the reconstruction techniques introduces 
significant artifacts in the intensity spectra. The speckle deconvolution (MOMFBD) has a rms noise in Stokes V/I of 0.32 % (0.20 %). 
The deconvolved spectra thus require a high significance threshold of about 1.0 % to separate noise peaks from true signal. A com- 
parison to spectra with a significantly higher signal-to-noise (S/N) ratio and to spectra from a magneto-hydrodynamical simulation 
reveals that the Gottingen FPI can only detect about 30 % of the polarization signal present in quiet Sun areas. The distribution of NCP 
values for the speckle-deconvolved data matches that of observations with higher S/N better than MOMFBD, but shows seemingly 
artificially sharp boundaries and unexpected changes of the sign. 

Conclusions. For our imaging data, both MFBD and speckle reconstruction are equivalent, with a slightly better and more stable 
performance of MFBD. For the spectropolarimetric data, the higher intensity contrast of the speckle deconvolution is balanced by 
the smaller amplification of the noise level in the MOMFBD at a comparable spatial resolution. The noise level prevents the de- 
tection of weak and diffuse magnetic fields. Future efforts should be directed to improve the S/N of the Gottingen FPI spectra for 
spectropolarimetric observations to lower the final significance thresholds. 

Key words. Sun: photosphere - Sun: spectropolarimeters 



1. Introduction 

One of the primary objectives of modern solar physics is the 
determination of the thermodynamic and magnetic structure of 
the solar atmosphere at the smallest spatial scales, including iso- 
lated magnetic elements in the quiet Sun and the internal fine- 
structure of larger features such as sunspots. An accurate deriva- 
tion of the atmospheric properties requires both high spatial and 
spectral resolution. Telescopes with apertu res of 1-1.5 m, suc h 
as the Swedish Solar Tel escope (SST, Scharm er et al] l20Q3ah . 
GREGOR (IVolkmer et al] 120071: b althasar et al] 12007b . or the 



New Solar Telescope (JDenker et a l. 2006), allow one to study 
the dynamics of the solar fine- structure at spatial scales down to 
50-100 km on the surface of the Sun. The next-generation solar 
telescopes of the 4-m class, i.e ., the Advanced Technology Solar 
Telescope (| Wagner et al. 112008k or the European Solar Telescope 
(Collados 2008), then will resolve the fundamental scale of the 
solar photosphere, i.e., the photon free path length. 

All observations from the ground are, however, degraded 
by the wavefront distortions caused by fluctuations of the re- 



fractive index in the Earth's atmosphere. To overcome this lim- 
itation, adaptive optics (AO) systems are used to correct the 
wavefront errors, but their compensation is only partial and the 
diffraction limit of the respective telescope is hard to reach. 
AO sy stems have been installed at the Dunn S olar Telescope 
(DST. iRimmele & Radickl [l998: Rimmele 200J), the German 
Vacuum Tower Telescop e (VTT , lvon der Luhe et alj|2003l) . and 
the SST (IScharmer ■ etd]l2003H). AO and multi- c onjugated AO 
system s JBeckersI 1 1981 [Berkefeld et al.1 120011: iMoretto et ail 
120041: lEllerbroek & Vog el 2009) with a larger corrected field of 
view (FOV) are also foreseen for all next-generation telescopes. 

For spectroscopic or spectropolarimetric observations taken 
with slit- spectrograph systems, it is difficult to improve the 
spatial resolution beyond the limit provided by the real- 
time correction of an AO system because of the essen- 
tially one-dimensional information at one position of the 
slit (see, e.g., [Ke ller & Johannesso nl 1 19951: Sutterlin ~& Wiehrl 
l200d:lAschwandenll2010l:lBeck et alJl201lL submitted to A&A). 
Contrary to that, instruments providing instantly a two- 
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dimensional FOV offer the possibility for a post-facto re- 
construction of the observations. Two-dimensional (2D) spec- 
trometers are commonly realized in the form of Fabry-Perot 
Interferometers (FPIs) that have a narrow tunable spectral band- 
pass and a high transmission. 

A first example of a 2D FPI spectrometer with three 
etalons w as designed for th e Au stralian solar faci li ty in 
Narrabri dRamsav et al.1 1 19701: lLoughhead et al.l 119781: iBravl 
Il988l) . Nowadays, almost all ground-based solar telescopes 
are equipped wi th FPI instruments: the Italian Panora mic 
Monochromator dBonaccini et al.l 1 19891: ICavallinil 1 19981) at 
the THEMIS telesco pe, the Interf erometric BIdimensional 
Spectrometer (IBIS ICavallini et al.l l2QQQfc ICavallinil [2006; 
iReardon & Cavallinil 120081) at theDS T, the visible-light and 
near-infrared imaging magnetographs (Den ker et alJ 120031) at 
the Big Bear Solar Observato ry, the CRisp Imaging Spectro 
Polar i meter (CRISP, e .g., Ivan Noort & Rouppe van der Voortl 
120081: IScharmer et all 120081) at the SST, and the Triple 
Etalon SQlar Spectrometer (TESOS, [K entisch eTet all 119981: 
Ivon der Ltihe & Kentischen l2QQ0fc iTritschler et alJ 120021) at the 
VTT. For the VTT, another 2D instrument was developed prior 
to TESOS by the Institu te for Astrophysics Gottingen, the 
so-called Gottingen FPI dBendlin et al.1 1 19921: iBendlinl 1 19931: 
iBendlin & Volkmedl 19931 1 19951) . This instrument will in the fu- 
ture be one of the post- focus instruments at the GREGOR tele- 
scope. Therefore, the spectrometer was largely renewed during 
the first half of 2005 (IPuschmann et al.ll2006|). Its de sign as the 
new "Gregor FPI" i s described in iPuschmann et al.1 (|2007) and 
iDenkeret all (120101) . 

Even if most of the 2D instruments were initially designed as 
spectrometers only, they now have all been upgra ded for (vector) 
spec tropolarimetry (see e .g., iBello Gonzalez & Kneer (2008) 
andlBalthasar et al.1 (l2009h for the Gottingen FPLlViticchieet all 
(2009b and lJudge et al.1 (l20Tol) for IBIS, and Bec k et al.l (12010) 
for TESOS). All 2D instruments use at least two channels, where 
one beam passes through the tunable FPIs with their high spec- 
tral resolution (narrow-band (NB) channel). The second channel 
uses only a broad-band (BB) interference filter to provide strictly 
simultaneous images with a high signal-to-noise (S/N) ratio for 
an accurate alignment of the individual wavelength steps and the 
application of post-facto reconstruction techniques. 

Two main reconstruction techniques have been developed 
and successfully applied to observation data up to now: the 
speckle reconstruction technique with the possibility of a subse- 
quent NB deconvolution of spectroscopic or spectropolarimetric 
data, and the blind deconvolution (BD) technique that can be ap- 
plied to both imaging and 2D spectropolarimetric data. 

For a speckle reconstruction, a fast series of about 50-100 
short-exposed images 7 bs(x, U) is taken at times t[. The observed 
images represent the image of the real object 7 r eai(x) modu- 
lated by the instantaneous optical transfer function (OTF) of the 
Earth's atmosphere and the telescope. In Fourier space, the cor- 
responding equation reads 



/ obs ( S ,^ = /real(s)-OTF(s,^), 



(1) 



where ~ denotes the Fourier transform, and x and s are the spatial 
coordinates and their corresponding Fourier equivalent, respec- 
tively. 

To recover the information of the real image, the Fourier 
amplitude and phases of / rea i are derived separately. The 
Fourier a mplitudes are usually determined by t he method of 
Labeyrie (1970) and the spectral ratio method (von der Lime 



19841) . Two main approaches exist for the estimate of 



masking algorithm (iKnox & Thompson! [19741: IWeigeltl Fl977l) . 
The large number of input images is necessary to ob- 
tain good atmospheric statistics. For the application to so- 
lar observations with low-contrast images, several specific 
spec kle reconstruct i on methods w e re developed (see, e.g., 
I von der Liihel 1 19931: Ide Boerl 1 19931: iMikurda & von der Luhd 
2006). Examples of recent speckle codes are the Kiepenheuer- 
Institute Speckle Interferometry Package (KISIP, Woger et al. 
2008]; IWoger & von der Luhe| |2008) and the Gottingen speckle 
reconstruction cod e (de Boer 1993). The latter was improved by 
IPusc hmann & Sailer (2006) to accoun t for the field - depen dent 
real-time correction by AO (see also iDenker et al.1 (120071) . or 
IWogerl (120101) for the KISIP code) a nd successfully applied 
in, e.g ., Puschmann & Wiehr (2006), Blanco Rodrig uez et al.1 
(120071) . or [Sanchez - Andrade Nuno et al.1 (120081) . For 2D spec- 
tro(polari)meters, the BB images can be directly speckle recon- 
structed because a sufficient number of images is available. The 
simultaneously recorded NB images - always only a few images 
per wavele ngth position - can then be deconvolved using the 
method of Kell er & von der Luhel (Il992l) estimating the instan- 
taneous OTF from the observed and rec onstructe d BB images 
(see a lso iKrieg et al.lfl999l: lJanssenll200l IBello Gonzalez et al.1 
120051) . 

In the case of the BD technique, one estimates both the real 
image and the OTF at the same time without using an OTF 
model based on the seeing statistics. This leads to a problem 
similar to the disentanglement of magnetic field strength and 
area filling factor for spectral lines in the weak-field limit: an 
infinite number of possible combinations can yield nearly the 
same result. The ambiguity is lessened by the use of multiple 
image frames, if the real image is assumed to be always the 
same during the observed fast image sequence. Compared with 
the speckle reconstruction, a significantly smaller number of im- 
ages, about 5-20, is sufficient (see, e.g.. Lofdahl et al .1120071) . 
Such a multi-frame blind deconvolution (MFBD, Schulz 1993; 
Ivan Kampen & Paxmanll 19981: lLofdahll 12002) can be applied to 
series of i mages of a single re al object. The approach was ex- 
tended bv lLofdahll (I2002J) and Ivan Noort et al.1 (12005b to cover 
also the case of multiple objects, e.g., simultaneous images of 
the same object in slightly different wavelengths or different po- 
larization states. This allows the reconstruction of the simulta- 
neously recorded BB and spectro (polari)metric NB images of, 
for instance, FPI instruments. The correspond ing multi-object 
multi-frame blind deconvolution (MOMFBD, Ivan Noort et al.1 
120051) code is freely available^. A drawback of the (MO)MFBD 
is the necessary computing effort that usually requires the use 
of a dedicated computer cluster or similar installations, whereas 
speckle codes can also be run within limitations on a single desk- 
top P C or nearly in real-time on a computer cluster (De nker et al.1 
l200ll) . 

Even if the application of reconstruction methods to imag- 
ing and 2D spectroscopic or spectropolarimetric data is a rou- 
tine task these days, there are only relatively few investiga- 
tions on the influence of the reconstruction on the data proper- 
tie s, e.g., [M ikurda et al. (2006) f or sp e ctrosc opic observations 
or Ivan Noort & Rouppe van der Voortl (120081) for MOMFBD- 
reconstructed data. The reconstruction introduces more severe 
changes of the data than, e.g., a simple flat fielding, and there- 
fore also may introduce systematic artifacts. 

FPI instruments sequentially scan the wavelength points; the 
independent observation and reconstruction of individual wave- 
length steps thus may compromise the spectral quality. This is 



the Fourier phases, i.e., the Knox-Thompson and speckle 
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Fig. 1. Setup of the observations in a top view. BS1: dichroic, 
splitting into NIR (> 800 nm) and VIS light. BS2: dichroic with 
its edge at 450 nm, splitting between the Gottingen FPI and the 
imaging channels. BS3: achromatic, 50/50. The two imaging 
channels in Can H and blue continuum used BB prefilters in 
front of the respective CCD cameras. 



even worse for spectropolarimetric observations, where addi- 
tionally the information on the polarization signal is encoded in 
the intensity of various modulation states. An independent scal- 
ing of the intensity in the individual images - varying across 
the FOV as well - therefore may easily introduce spurious po- 
larization signals because they are finally derived from intensity 
differences. 

In this contribution, we investigate the performance of 
speckle and BD reconstruction techniques. Because our observa- 
tions provided us with both imaging and 2D spectropolarimetric 
data, we extended the study to both types of data. In Sect. [21 we 
describe the observational setup and the data acquisition. Section 
[3] explains the data reduction and reconstruction methods for the 
individual data sets. The results presented in Sect. |4] are summa- 
rized and discussed in Sect. Section [6] provides the conclu- 
sions. 



2. Observations 

The observations consist of a time-series of 65 minutes taken 
on 2009 June 7 between 07:45 and 08:50 UT at the VTT with 
real-time correction of seeing distortions by the K iepenheuer- 
Institute AO system (KAOS, Ivon der Ltihe et alJ [2003). The 
Fried parameter r$ as estimated by the KAOS system was on av- 
erage about 10 cm during the observations. The target area was 
a quiet Sun (QS) region at disk center. 



2.1. Setup 

We used the setup prepared for the ground-based support ob- 
servations during the Sunrise flight (see Fig. [B, following the 
co ncept for mu lti- wavelength observations at the VTT described 
in iBeck et aD (I2007I) . The light from the telescope was first 
split into its visible (VIS) and near-infrared (NIR) parts by a 
dichroic beamsplitter (BS1) with a separation edge at 800 nm. 
The NIR light (> 800 nm) was transmitted towards the Echelle 
spectrograph, where the Tenerife Infrared Po larimeter (TIP, 
Mar tinez Pillet et alJ[l999l: ICollados et alJl2007h was mounted. 
TIP took high cadence spectra of the 1083 nm range at a fixed 
slit position. These data, however, will not be used in the current 
study. The VIS part was reflected towards the Gottingen FPI. 
Two imaging channels in blue continuum (BC) and Ca n H were 
fed from the VIS beam by a dichroic pentaprism (BS2) that re- 
flects all light below 450 nm to the side. The reflected light was 
evenly shared between the two imaging channels using an achro- 
matic 50/50 beamsplitter (BS3). For the final wavelength selec- 



tion, BB interference filters (431.3+0.3 nm and 396.85+0.3 nm) 
were mounted directly in front of the two imaging cameras. 

2.2. Data properties 

Imaging data The two imaging channels used PCO 4000 cam- 
eras in a 2x2 binning mode yielding a spatial sampling of 07071 
per pixel for the blue continuum channel and 07076 per pixel 
for Can H. The FOV was about 110" x 80". The cameras were 
run independently of the rest of the instruments in a streaming 
mode. Because the disk write time was different due to the im- 
age size and the intrinsic write time of the PC hard disks, we 
added a short delay time to the faster channel to obtain nearly 
identical frame rates. We then finally obtained bursts of 100 (96) 
images each 17 seconds with an exposure time of 7.5 (20) ms 
for the blue continuum channel (Ca n H). During the first 50 min 
of uninterrupted observations, we collected 200 bursts in each 
channel that were individually reconstructed. 

Gottingen FPI data The vector-polarimetric mode of the 
Gottin gen FPI in its l atest setup at the VTT was used 
dPuschmann et al.1 120061 : iBello Gonzalez & Kneed l2008h . We 
scanned the Fei line at 630.25 nm in 28 steps of 2.31 pm step 
width in wavelength, covering a range from about 630.23 nm 
to 630.29 nm. At each wavelength step, four modulation states 
were acquired with an exposure time of 20 ms for each image. 
The modulation sequence was then repeated eight times to have 
a sufficient number of images for the later deconvolution. Strictly 
simultaneous BB images (28x4x8 = 896 images per wavelength 
scan) were taken on a second camera with a BB prefilter cen- 
tered at 630 nm. We chose a frame rate of 10 Hz for the ob- 
servations because all higher frame rates led sporadically to a 
complete loss of the synchronization between the modulation by 
the liquid crystals and the camera read-out. The cadence for a 
line scan was therefore limited to 90 seconds. The spatial sam- 
pling was 0711 per pixel across a FOV of 38" x 58" before the 
reconstruction. The final size of the FOV then slightly varies de- 
pending on the reconstruction method. 

2.3. Additional data 

For comparison with the spectra from the Gottingen FPI, we 
will use three additional data sets. One is an observation t aken 
with the Hinode spectropolarimeter (S RlKosugi et al.ll2007l) that 
is described in more detail in Pus chmann et al.l (l2010alfbh . This 
scan with the Hinode/SP was taken on the I st of May 2007 
and used an integration time of 4.8 s per scan step of about 
0715 width, covering an active region close to disk center. The 
root-mean- square (rms) noise level in Stokes V in a continuum 
wavelength window was about 1.2 x 10" 3 of I c . We cut out 
a 20" x 20" QS section of its FOV for this study. The sec- 
ond observation used was taken with th e POlarimetric LIttrow 
Spectrograph (TOLIS. lBeck et al.ll2005bl) . We selected an obser- 
vation of QS on disk center obtained on the 29^ of Aug 2009, 
cutting out again a 20" x 20" section. The integration time per 
scan step of 075 width was 26 sees, yielding a rms noise of about 
3 x 10" 4 of I c (see the POLIS data archivqj for more details). 
These data provide the best reference for the polarization signals 
to be expected in QS because of their low noise level. As third 
data set, we used spectra derived from a magneto-hydrodynamic 



2 http://www3.kis.uni-freiburg.de/~cbeck/POLIS_archive/ 
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(MHD) simulation. T he MHD simulation was performed with 
the C0 5 BOLD code (Frevtag et al. 2002) and is described in 
more detail in Schaffenberger et al. (2005, 2006). For one simu- 
lation snap shot, spectra of the two Fe i lines at 630 nm were syn- 
thesized with a spectral sampling of 2 pm per pixel (for m ore de- 
tails see lSteiner et al.ll2008l: [Beck & Rammachedl2010l) . Figure 
[T2l later on shows maps of the continuum intensity I c for the 
three additional data sets. 



3. Data reduction 

3.1. Imaging data 

The data from the two imaging channels (and the Gottingen 
FPI) were corrected for flat field inhomogeneities and the dark 
current. For the imaging data, we then applied two reconstruc- 
tion methods to every burst of about 100 images: the mod ified 
versio n of the "Gottingen" speckle code dPuschmann & Sa iler 
2006) and the MFBD met hod. For the latter, w e used the ver- 
sion of the algorithm of Ivan Noort et al.l (120051) installed at the 
Instituto de Astrofisica de Canarias (IAC) that is prepared for 
parallel computing. 



3.2. Gottingen FPI data 

The spectra of the Gottingen FPI consist of a set of NB inten- 
sity images I(A(, rrij, k), where the index / = 1, . . . , 28 denotes 
the number of wavelength steps, rrij the four modulation states 
j = 1, . . . , 4, and k = 1, . . . , 8 the index of repetitions of the 
modulation cycle at one wavelength step. A simultaneous set 
of BB images B(At, rrij, k) is available from the second camera 
of the Gottingen FPI. The two orthogonally polarized beams in 
the NB images (I + and /") and the corresponding FOV of the 
BB images were extracted and aligned with pixel accuracy to 
avoid an interpolation of the data prior to any image reconstruc- 
tion. The Gottingen FPI spectra were then reduced with three 
different methods: simple destretching, speckle deconvolution, 
and MOMFBD. 



Destretching For destretching and later also the speckle decon- 
volution, we first speckle-reconstructed the 896 BB images of 
each wavelength scan to obtain a reference image # spe ckie- For 
destretching, we then determined the necessary shifts in the FOV 
to co-register the individual BB images during the wavelength 
scan with the speckle reference image and applied the same cor- 
rections to the NB images P (At, rrij, k). We then averaged the 
co-registered eight repetitions of each modulation state to get 
/*(/!/, rrij) =1/8 ZLi ^(^ m h ^)- Destretching and speckle de- 
convolution were applied separately to the two beams I + and 
/", whereas the MOMFBD deconvolves both beams at the same 
time (see Eq. 0]) below). 



Speckle deconvolution For the reconstruction of the spectropo- 
larimetric data we adapted the Got tingen NB deconvolution code 
dKrieg et al.l 1 19991: iJanssenl [2003) for an application to vector 
polarimetric measurements. The eight individual NB images 
I(Af, rrij, k) of each modulation state rrij and each wavelength \ t 
are deconvolved using the information obtained from the simul- 
taneous BB images and the speckle-reconstructed reference im- 
age ^speckle- In the Fourier domain, the corresponding equation 



reads 



f±n . „ Xl i n*i,m J9 k)B*(A i9 m J9 k) 

I (A h rrij) = F — £ speck i e 

Y? k=l \B{A u mjM 2 



(2) 



where ~ denotes the Fourier transform, * t he complex conjugate, 
and_F is a noise filter (optimum filter, see Kell er & von der Luha 
Il992l) . The inverse Fourier transform then yields the wanted 
/*(/[;, rrij). The deconvolution is applied individually to overlap- 
ping subfields of 64x64 pixels because Eq. ([2] ) is only valid in- 
side an isoplanatic patch (Ivon der Luhdl 1993b . 



MOMFBD For the MOMFBD, all images from a single wave- 
length scan were fed simultaneously into the code to obtain one 
reconstructed BB image # rea i and one image per modulation 
state, wavelength position, and orthogonally polarized beam, 
PiAi.mj). We used the default settings of the MOMFBD code 
which employs a decomposition of the wavefront phases into 
36 Karhunen-Loeve modes and an equal weight for NB and BB 
data. The code then minimizes 

28 4 8 

L = Z Z Z Z |; ^' m > k) - £* ( *' m ^ • OTF(/ ' j' * )|2(s) 

;=0 7=0 k=\ S 

+|/- bs (^-, rrij, k) - /- al (^-, rrij) • OTF(/, j, k)\\s) 



+ |fi obs Wi, m P k ) ~ ^real ' OTF(l, j, k)\\s) , 



(3) 



where s denotes the Fourier equivalent o f the spatial coordi- 
nates (lLofdahlll2002l:lvan Noort et al.ll2005h . This reconstruction 
is also performed for subfields of 64x64 pixels. 

As became clear during the subsequent investigation of the 
deconvolved spectra, the solution with equal weight for the BB 
and NB channel in the MOMFBD is only one possible option 
with specific characteristics. We run a series of additional test 
deconvolutions that are described in Appendix |Al but finally de- 
cided to stick to the initial settings because they provide the low- 
est noise level in Stokes V at continuum wavelengt hs of all test 
runs and maintain the spectral quality best (see Sect. IA.5l) . 

All three reduction methods for Gottingen FPI spectra yield 
finally a set of 28x4 intensity images for the two beams, 
/*(/!/, rrij). The spectra were then corrected for the transmission 
curve of the NB prefilter, using a separate wavelength scan ob- 
tained with a continuum light source. Afterwards, we applied a 
correction for the field-dependent blue shift in the spectra caused 
by the collimated mounting of the etalons. The necessary blue 
shift correction was determined from the line-core position of 
the solar Fei line in the average flat field frames. The polar- 
ization information could now be retrieved by the polarimetric 
calibration procedure described in the next section. 

3.3. Polarimetric calibration 

The polarimetric calibration of the Gottingen FPI data was car- 
ried out in a two-step procedure with a correc tion for the instru- 
mental polarization induced by the telescope (JBeck et al.l l2005a) 
and the cross-talk caused by the optics behind the instrumen- 
tal calibration unit (ICU) of the VTT. The calibration data ob- 
tained with the ICU were used to determine the polarimeter re- 
sponse function with the gen eri c approach for the VTT described 
in, e.g.. lBecket~aTl(l2005bl) or lBecketail (l2010h . The efficien- 
cies for measuring the components of the Stokes vector (cf. 
del Toro Iniesta & Collados 2000) were about (e/, €q, eyey) = 
(0.93,0.56,0.49,0.42) in our case, with a total efficiency for 
measuring polarization of e po i = 0.85. 



K.G.Puschmann & C.Beck: Application of speckle and MOMFBD/MFBD techniques on 2D data 



unfiltered 



low-pass filtered 





Table 1. Rms noise in Stokes V/I in percent. 
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unfiltered filtered running mean 



v / i [%] 




0.0 0.2 0.4 0.6 0.8 1.0 

^Nyquist 



Fig. 2. Histograms of the polarization amplitude of Stokes V/I 
in a continuum window. Top left: unfiltered data for speckle de- 
convolution (black), MOMFBD (red), and destretching (blue). 
Top right and bottom left: the same for application of a low- 
pass filter and of a weighted running mean, respectively. Bottom 
right: Fourier power vs. spectral frequency. Solid/dotted/dashed 
lines: unfiltered/low-pass filtered/running mean. The vertical 
black line corresponds to a wavelength of 10 pm. 



Residual / — » QUV cross-talk was removed by forcing the 
continuum polarization level to zero. We determined the average 
ratio a quv = QUV /I in the continuum wavelength range from 
630.275 to 630.291 nm separately in each profile and then sub- 
tracted a quv ' IW from Stokes QUV(A). This approach com- 
pensates automatically any possible variation of the cross-talk 
level across the FOV, but could only be used because contin- 
uum wavelength points without intrinsic polarization signal were 
available. Scanning only a smaller wavelength range around the 
solar spectral line does not permit to use this correction because 
it cannot be assumed that the polarization signal at these wave- 
lengths should always be zero in real observations. 

3.4. Spectral noise filtering 

The reconstruction process changes the intensities at different 
spatial locations and at different wavelength points indepen- 
dently of each other. Because the polarization signal is derived 
from intensity differences between the four images belonging 
to the different modulation states, a relative change of intensity 
by the reconstruction process can produce spurious polarization 
signals. Imperfections in the spatial alignment of different mod- 
ulation states can also have the same effect. The spectra of the 
speckle-deconvolved data show a large noise level in the polar- 
ization signal of the Stokes V/I signal at continuum wavelengths 
(630.275 nm to 630.291 nm) with peak values larger than 1.5 % 
(top left panel of Fig. [2]), which requires some sort of filtering. 
The noise level in the MOMFBD or destretched data is lower 
than that of the speckle-deconvolved data by a factor of about 
1.5; the MOMFBD data show a slightly higher noise level than 
the destretched data (first column of Table [T]). 

We used two different approa ches for the spectra l noise 
filtering: a low-pass filter as in Bello Gonza lez et al. I (2009) 



Speckle deconvolution 0.50 0.32 0.32 

MOMFBD 0.32 0.20 0.20 

Destretching 0.29 0.18 0.18 



and a weighted running mean in the spectral dimension over 
three wavelength steps (= 6.9 pm). For the latter, we calcu- 
lated the new values of Stokes I QUV at the wavelength X\ by 
2 j jl + _\ a j+1 IQUV(A i+j ) with a= (0.25,0.5,0.25). The low-pass 
filter was set to cut all power for frequencies above half of the 
Nyquist frequency of (2 x 2.31 pm)" 1 . The upper right and lower 
left panels of Fig. |2] show the histograms of the Stokes V contin- 
uum polarization signal after applying the two filters to the spec- 
tra. In both cases, the peak values of the noise reduce to about 
1 % for the speckle deconvolution and to about 0.5% for the other 
two data reduction methods. 

In the lower right panel of Fig. [2 the Fourier power as a 
function of the spectral frequ ency normalized with the N yquist 
frequency is displayed (cf. iBello Gonzalez et all 120091 their 
Fig. 1). The power spectra show that the main power is concen- 
trated at small frequencies corresponding to a large wavelength 
extent, but the power in the unfiltered data stays at about 50% 
of the maximum value up to the frequency limit because of the 
high noise level. The weighted running mean has a very similar 
effect as the low-pass filter and reduces the power to about 20% 
at half of the Nyquist frequency. For comparison, the full width 
at half maximum (FWHM) of the Fe i line at 630.25 nm of about 
10 pm is marked with a vertical black line. Because any spectral 
features should be broadened to a comparable extent by thermal 
and Doppler broadening, none of the two filters should thus have 
removed significant solar information. 

Figure [3] shows example spectra of Stokes I QUV from a lo- 
cation inside the network (top panel) and a small isolated polar- 
ization patch (bottom panel) for all three data reduction meth- 
ods before and after the application of the two filters. The low- 
pass filter and the spectral running mean produce identical spec- 
tra. The largest noise peaks in the speckle-deconvolved data are 
in some cases not fully removed by the filtering, but spectrally 
broadened such that they resemble a false solar polarization sig- 
nal with an amplitude of up to 1 % of Stokes / in some cases 
(e.g., in the black rectangle in Q in the uppermost row of Fig. [3]). 
Some of the peaks, e.g., in Stokes V of the speckle-deconvolved 
data in the lower panel of Fig. [3] are located at wavelengths that 
would require Doppler shifts of several tens of kms -1 or field 
strengths far exceeding 2 kG. These features are also missing 
in the other two reduction methods. This excludes a solar ori- 
gin, but indicates that it is not straightforward to separate real 
from false polarization signals at the noise level of the speckle- 
deconvolved data (see Sect. l4.3l below). 

The shape of the intensity profile stays virtually the same for 
all data reduction methods and applied filters. Table [T] lists the 
rms noise in Stokes V/I in a continuum window before and after 
noise filtering. Because the weighted running mean and the low- 
pass filter have nearly the same effect on the spectra, we selected 
the weighted running mean as our method of choice for noise 
filtering. It should have the lowest impact on the spectral proper- 
ties because it only operates on directly neighboring wavelength 
steps. 



K.G.Puschmann & C.Beck: Application of speckle and MOMFBD/MFBD techniques on 2D data 



Speckle deconvolution 




unfiltered 

low — pass filt. 

running mean 



07 








0.01 




II JS 




0.00 


*yfy 


T 


V 'l 


0.01 




I 




1)7 










0.25 0.27 0.29 
A-630 [nm] 



0.25 0.27 0.29 
A-630 [nm] 



0.23 0.25 0.27 0.29 
A-630 [nm] 



0.23 0.25 0.27 0.29 
A-630 [nm] 




unfiltered 

low — pass filt. 

running mean 




0.02 
0.01 - 
0.00 
-0.01 




0.25 0.27 0.29 
A-630 [nm] 



Destretching 



0.23 0.25 0.27 0.29 
A-630 [nm] 



0.23 0.25 0.27 0.29 
A-630 [nm] 



0.23 0.25 0.27 0.29 
A-630 [nm] 




0.02 
0.01 - 
0.00 
-0.01 



0.05 
0.00 



0.25 0.27 0.29 
A-630 [nm] 



0.23 0.25 0.27 0.29 
A-630 [nm] 



0.23 0.25 0.27 0.29 
A-630 [nm] 



0.23 0.25 0.27 0.29 
A-630 [nm] 



Speckle deconvolution 




- unfiltered 

- low-pass filt. 

- running mean 



0.02 
0.01 










0.00 


st 






-0.01 
-0.02 












0.03 




0.02 




r\ 


0.01 








0.00 








0.01 








0.02 





V : 


0.03 





0.23 0.25 0.27 0.29 
A-630 [nm] 



0.23 0.25 0.27 0.29 
A-630 [nm] 



0.23 0.25 0.27 0.2 
A-630 [nm] 



0.23 0.25 0.27 0.29 
A-630 [nm] 




- unfiltered 

- low-pass filt. 

- running mean 



0.23 0.25 0.27 0.29 
A-630 [nm] 



0.23 0.25 0.27 0.29 
A-630 [nm] 



Destretchin 


g 


unfiltered 

low-pass 


r 

filt. 

lean" 



0.02 

0.01 

0.00 

-0.01 

-0.02 
0. 

0.02 
0.01 
0.00 
-0.01 
-0.02 



0.03 
0.02 




0.23 0.25 0.27 0.2 
A-630 [nm] 



0.23 0.25 0.27 0.29 
A-630 [nm] 



0.03 




0.02 




0.01 


/^ A 


0.00 


^-"V^H 


0.01 




0.02 




0.03 





0.23 0.25 0.27 0.29 
A-630 [nm] 



0.23 0.25 0.27 0.29 
A-630 [nm] 



0.23 0.25 0.27 0.2 
A-630 [nm] 



0.23 0.25 0.27 0.29 
A-630 [nm] 



Fig. 3. Example spectra in the network (top 
panel) and in a small isolated polariza- 
tion patch (bottom panel). Left to right: 
Stokes IQUV. Top to bottom: speckle de- 
convolution, MOMFBD, and destretched data. 
Black/red/blue lines: unfiltered data/low-pass 
filter/running mean. The black rectangles de- 
note noise peaks that could be mistaken for so- 
lar polarization signal after the filtering. 



4. Results 

4.1. Overview maps 

For a first visual inspection of the results, Figs. |4] and [5] show 
intensity maps of the reconstructed imaging data (blue contin- 
uum, Ca ii H, 630 nm BB data) and all Stokes parameters of the 
spectropolarimetric data, respectively. 



Imaging data In case of the imaging data in the blue contin- 
uum and Ca (upper two rows in Fig. 0]), the two reconstruction 
methods yield generally similar results with only a slightly dif- 
ferent behavior in the two wavelength channels. For this specific 
burst, the MFBD of the Ca data has a slightly better spatial res- 
olution as can be seen when comparing the shape of the bright- 
enings around (x,y) « (15", 8") in the speckle and MFBD data 
that is displayed magnified on the right-hand side of Fig.ffl The 
MFBD resolves three individual round brightenings inside the 
red square in the magnification, whereas the speckle reconstruc- 
tion shows the "blooming" effect, i.e., a radial smearing of iso- 
lated bright features. This effect produces a bright ring around 



a central darkening in the difference image. Some more exam- 
ples throughout the FOV are marked by the small grey squares in 
the Ca difference image. The burst shown was randomly picked 
near the end of the time- series, but a comparison of the other 
reconstructed Ca images showed always the same trend, i.e., a 
slightly better performance of MFBD in Ca in terms of spatial 
resolution. For the blue continuum data (middle row), the re- 
construction results of both methods are virtually identical. The 
structures in the difference image are usually located at the edges 
of granules and come from the slightly different placing of indi- 
vidual subfields on a subpixel scale. For the reconstruction of 
the BB data at 630 nm taken with the Gottingen FPI (bottom row 
of Fig. [4]), the results differ unexpectedly strong. The speckle re- 
construction shows a significantly larger intensity contrast than 
the MOMFBD. The contrast difference reproduces the complete 
granulation pattern in the difference image (see Fig. lA.2l for the 
trade-off between noise level and contrast in the MOMFBD). 



Spectropolarimetric data For the spectropolarimetric data of 
the Gottingen FPI, Fig. [5] shows maps of the four Stokes pa- 
rameters. The map of the continuum intensity I c was taken 
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Fig. 4. Overview maps of the imaging data. Bottom to top: Gottingen FPI BB channel, blue continuum, and Ca. Left to right: 
speckle reconstruction, (MO)MFBD, and difference image. The display ranges for the intensity maps are 0.8 < I < 1.5 for Ca 
and 0.8 < / < 1.2 otherwise. The difference images are thresholded at ±0.05. The large white rectangles denote the region shown 
magnified in the most right column, the small grey squares indicate locations with ring structures in the Ca difference image. 



at the first wavelength step of the spectral scan. For Stokes 
QUV, the maps of the wavelength-averaged unsigned polariza- 
tion signal, f \QUV/I\(A)dA/ JdA, are displayed. As for the BB 
data, the speckle deconvolution has a significantly larger con- 
trast than MOMFBD in I c . No location inside the observed 
FOV shows a significant linear polarization signal (second and 
third row of Fig. [5]). The noise level in the polarization states 
as given by the background patterns in QUV is smallest for 
the destretched data set and largest for the speckle-deconvolved 
data. The latter also show a weak grid pattern at the noise 
le vel with an edge length of ab out 3", which was also seen 
in lSanchez-Andrade Nun o (2009). The reduced intensity of the 
grid pattern in the wavelength- averaged polarization signal im- 
plies a reduction of the noise le vel on those pixels where indi - 
vidual subfields were connected ( Sanchez- Andrade N uno 2009). 
The size of the grid matches the settings of the speckle recon- 
struction, i.e., subfields of 64 pixels (= 7.04") with a 33 % over- 
lap to each neighboring one, where the apodization effects a 
gradual change from one subfield to the next. 



Table 2. Continuum contrast in percent. 



The polarization signals in Stokes V show some differences 
between the different data reduction approaches. Several small- 
scale (roundish) features are well defined and resolved from each 
other in the speckle-deconvolved and MOMFBD data (e.g., in- 
side the two white rectangles, especially in the lower part of 
rectangle b), whereas they appear significantly smeared out in 
the destretched data. The destretched data, however, show a po- 
larization signal at the same locations as the deconvolved data 
sets. This implies that both deconvolution approaches in general 
do not reveal "new" signal that is absent in the destretched spec- 
tra, but only localize the origin of the polarization signal more 
precisely by removing the spatial smearing. 



Type 


BB 


NB 


BC 


CA 




630 nm 


630 nm 


431 nm 


396 nm 


Speckle 


8.0 


6.3 


7.6 


6.4 


(MO)MFBD 


4.6 


4.3 


7.8 


6.5 


Destretching 


- 


3.6 


- 


- 


MHD 


- 


15.3 


- 


- 


SP 


- 


7.1 


- 


- 



4.2. Continuum contrast and spatial power spectra 

We estimated the spatial resolution for the various maps quan- 
titatively using the rms contrast and the power spectra of the 
intensity maps. The contrast gives only a partial estimate of the 
spatial resolution because it depends on various influences. Table 
[21 lists the rms contrast values. The numbers support the visual 
impression of the previous section: for the imaging channels, 
MFBD and speckle reconstruction give nearly identical con- 
trasts. For the continuum intensity of the spectra or the BB chan- 
nel, the contrast of the speckle deconvolution or reconstruction is 
larger by 2-4% in absolute numbers (see Appendix lAl for other 
MOMFBD solutions). The observations have not been corrected 
for stray light yet, which will increase the absolute contrast val- 
ues but presumably not their differences. Table |2] also lists the 
corresponding continuum intensity contrast in the MHD simula- 
tion and the Hinode/SP data (see Sect. 12.3b for comparison. 

The spatial power spectra are shown in Fig.O They are nor- 
malized individually to the power at the first non-zero frequency. 
The imaging data in blue continuum and Ca show a similar be- 
havior, i.e., the power of the speckle reconstruction drops faster 
than the MFBD at spatial scales of about 075 and levels off to a 
constant level earlier (see also Matson 2008). This trend is more 
pronounced for the Ca data, as expected from the visual impres- 
sion (Fig. [4]). For the spectra of the NB channel, the difference 
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Fig. 5. Overview maps of the NB data. Left to right, speckle deconvolution, MOMFBD, and destretching. Bottom to top: 0.85 
< I c <1.15, wavelength-averaged unsigned Stokes Q < 2%, U < 2%, and V < 2%. The polarization maps are displayed on a 
logarithmic scale. 




Fig. 6. Spatial power spectra. Thick red: (MO)MFBD data. Thin 
black: speckle-reconstructed/deconvolved data. Top to bottom: 
blue continuum, Ca n H imaging data, continuum window in the 
630 nm spectra, 630 nm BB data (each vertically displaced by 
two units). The thick blue triple -dot dashed, dash- dotted purple, 
and dashed green lines show the power spectra of Hinode/SP 
630 nm spectra, a MHD simulation, and the destretched spectra, 
respectively. 



between the speckle deconvolution and the MOMFBD is small. 
The speckle deconvolution shows slightly larger power at spa- 
tial scales of 1" and for scales smaller than 074, where the lat- 
ter produces the higher noise level of the speckle deconvolution. 
The power spectrum of the destretched spectra (dashed green 
line) lies below that of both deconvolved spectra, implying both 
a lower noise level and a worse spatial resolution. Both the MHD 
simulation (red dash-dotted) and the data from the Hinode/SP 
(thick blue) show a power spectrum with a nearly linear shape in 
the logarithmic plot that the power of the speckle-reconstructed 
BB images matches for spatial scales larger than about 075. 



The effective spatial resolution is usually determined from 
the location where the power curve changes its slope signifi- 
cantly and levels off to a roughly constant level. The exact loca- 
tion, however, is difficult to pinpoint in observed power spectra. 
We would attribute a spatial resolution of about 0725 (= v = 4) 
to the imaging data and about 074 to the deconvolved NB spec- 
tra, whereas the destretched NB spectra reach rather only about 
076. The power spectrum of the Hinode/SP data shows no real 
change of the slope, implying a diffraction-limited performance 
with 0732 as the theoretical resolution limit. 
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Fig. 7. Average profiles of the Gottingen FPL Black line 
with crosses', speckle deconvolution. Red line with asterisks: 
MOMFBD. Blue line: destretched data. The thick green (thin 
green line with diamonds) line shows the corresponding FTS at- 
las profile before (after) the convolution with the instrumental 
profile. 



4.3. Spectroscopic properties 

To investigate a possible effect of the deconvolution methods on 
the shape of the spectra, we determined the following parame- 
ters: the average profile, the line width, the relative line asymme- 
try, the line-core intensity, and the line-core position (converted 
to the corresponding LOS velocity). 

Averaging over the FOV yields nearly identical profiles for 
all data reduction methods (Fig. [7]). Because the deconvolution 
should preserve the total intensity, and therefore also the av- 
erage profiles, the very small differences imply that all three 
data reduction methods performed correctly in this respect. 
We matched the observ ed spectra and the corresponding sec- 
tion of the FTS atlas (Kuruczet al.|ll984|) with the approach 
used i n lAUende Prieto et al.l (12004) and lCabrera Solatia et al.l 
(120071) using a stray light offset a and a subsequent convolu- 
tion of the FTS spectrum with a Gaussian of width cr. The 
Gaussian curve is assumed to account for the spectral resolution 
of the instrument. The best-fit was reached for (a,o~) = (14 %, 
1.65 pm), which yields a FWHM of the Gaussian of 3.88pm 
and thus a spectral resolution of about 162.000, somewhat be- 
low the theoretical li mit of the Gottingen FPI of about 240.000 
dBello Gonzalez & Kneerfl2008l) . 

Prior to the determination of the other line parameters, we 
removed the telluric line blend at 630.27 nm from the spectra by 
fitting a Gaussian to it and subtracting the corresponding contri- 
bution from the profile. The line width was provided by a sub- 
sequent Gaussian fit to the solar line. The line-core position of 
the solar line was derived using a Fourier method for the close 
surroundings of the intensity minimum. The position was finally 
converted to the corresponding LOS velocity using the average 
line-core position as zero reference. The line asymmetry was de- 
termined as the difference of the area between the continuum in- 
tensity level and the line profile to the red and the blue from the 
core position for an equal distance of ± 16.2 pm in wavelength. 
The difference was then normalized with the total area. Prior 
to the determination of the line asymmetry, the whole profile 
was shifted to have the line-core position exactly at one of the 
observed wavelength steps to avoid spurious asymmetry values 
caused by the discrete wavelength sampling. The line-core in- 
tensity was derived from the minimum intensity in the spectral 
line. 
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Fig. 8. Line parameters. Top to bottom: line width, line asymme- 
try, line-core intensity, and line-core velocity. Left to right: de- 
stretched data, speckle-deconvolved data, MOMFBD data. Tick 
marks are in arcsec. The display ranges are given by the gray- 
scale bar at bottom. 



Figure [8] shows maps of the retrieved line parameters. The 
image orientation is at 90 degrees to the previous figures. The 
line-core intensity and line-core velocity show a similar pat- 
tern for both deconvolution techniques and a considerably im- 
proved resolution when compared with the destretched data set. 



K.G.Puschmann & C.Beck: Application of speckle and MOMFBD/MFBD techniques on 2D data 



Table 3. Line parameters of the intensity profiles and their rms 
values. 



Data 


line width 


line asym. 


core int. 


Vlos rms 




pm 


% 


%ofI c 


kms -1 


Destretching 


12.33+0.71 


1.5+2.4 


46.1+2.2 


0.296 


MOMFBD 


12.47+0.57 


2.1 + 1.9 


45.7+2.6 


0.347 


Speckle deconv. 


12.57+0.69 


2.7+2.6 


45.5+3.2 


0.454 


Hinode/SP 


11.71+0.99 


-0.2+1.4 


38.8+5.1 


0.537 


MHD Sim. 


12.18+1.21 


-0.1+2.5 


38.3+7.6 


1.029 


MHD, degr. 


12.43+0.76 


-0.6+1.2 


40.2+5.0 


0.728 



Locations of strong magnetic fields can be clearly identified in 
the line-core intensity (see also Pusch mann et al.ll2007h and the 
line- width map. For both the line asymmetry and the line width, 
the deconvolution changed, however, to some extent the spatial 
structuring across the FOV when compared with the destretched 
data set. A weak grainy pattern can be identified in the line width 
of the speckle-deconvolved da ta that was noted as typ ical for this 
reconstruction method also by Mikurda et al. (2006). 

Table [3] lists the average values of the line parameters 
and their rms fluctuations for the Gottingen FPI spectra, the 
Hinode/SP data set, and the MHD simulation. The table is 
roughly sorted by increasing spatial resolution from top to bot- 
tom (cf. Sect. 14. 21) , except for the last line that corresponds to the 
MHD simulation after a spatial degradation with a 0752 wide 
Gaussian kernel, which intends to mimic the resolution of the 
POLIS data. For most quantities, e.g., LOS velocity, line width, 
and line-core intensity, a clear trend of increasing rms variations 
can be seen with the spatial resolution. The larger line width and 
higher line-core intensity of the Gottingen FPI spectra in com- 
parison to the other two data sets are presumably caused by the 
instrumental profile and the resulting lower spectral resolution 
of the Gottingen FPI. 

4.4. Spectropolarimetric properties 

The intensity spectra showed no strong effect by the deconvo- 
lution, but the spectra still have a considerable noise level even 
after the spectral filtering. We thus first focus on quantities that 
allow to determine the significance level of the polarization mea- 
surements like the Stokes V amplitude and the polarity of the 
magnetic fields before we turn to the net circular polarization 
(NCP). 

We used the maximum value of the unsigned polarization 
signal in Stokes V as proxy for the presence of a polarization 
signal because the solar Stokes Q and U signals are below the 
noise. The polarity was defined as the order of maximum and 
minimum signal in the signed Stokes V profile and set to ±1 
accordingly. We used a threshold of 2.5 times the rms noise level 
cr in the Stokes V spectra (see the last column of TableQ} as 
initial criterion for a significant signal (cf . the 2 cr level used in 
parts of iBello Gonzalez & Kn eer 2008). The polarity was set to 
zero for all locations with V max < 2.5 cr. Figure [9] shows the 
maps of the Stokes V amplitude, the polarity, and the mask of 
locations with V max > 2.5 cr (black shading). 

From a comparison of the masks or the polarity maps with 
the Stokes V amplitude, it is clear at a first glance that the 2.5 cr 
threshold is insufficient because in many places no significant 
corresponding signal can be seen in the Stokes V map. The po- 
larity of the destretched and MOMFBD data matches for most 
locations where the signal level in the V-amplitude map ex- 
ceeds the display threshold of 1.5 % (e.g., inside the black rect- 




Fig. 9. Left to right, maximal V amplitude thresholded at 1.5%, 
polarity, locations with V max > 2.5 cr (black), and polarity 
of the significant polarization signals with the final rejection 
thresholds. Top to bottom: destretched data, MOMFBD, speckle- 
deconvolved data. The white rectangle marks the area that is dis- 
played magnified in Fig.[l0l The two vertical dash-dotted lines 
denote the location of the two cuts shown in Fig.fTTl 
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Fig. 10. Magnification of the area marked in Fig.[9j Left to right. 
continuum intensity, maximal V amplitude < 1.5%, polarity for 
V > 2.5 cr, and polarity with the final rejection thresholds. Top: 
MOMFBD. Bottom: speckle deconvolution. 



angles). The speckle-deconvolved data sometimes show a po- 
larity opposite to the other two data sets (e.g., at some places 
inside the lower black rectangle). The polarities in both decon- 
volved data sets tend to fluctuate on small spatial scales at the 
borders of large-scale con nected polarization patches (cf. also 
Marti nez Pillet et alJl201 ll their Sect. 9.3). 



10 



K.G.Puschmann & C.Beck: Application of speckle and MOMFBD/MFBD techniques on 2D data 



I MOMFBD M)ESTR ^SPECKLE ^MOMFBD 




'polis Vpous + noise 


40 






; i » 




30 








C: , 


IB 


'°'m 




10 


ij 


i- ^ '% 


j_ 


, ,%^ -r> 





-1515 -1515 -1515 -1515 
pm pm pm pm 



-1515 -1515 -1515 -1515 
pm pm pm pm 



-1515 -1515 -1515 
pm pm pm 



Fig. 11. Stokes / and V spectra along the 
two cuts marked in Fig. [9] (left two pan- 
els), and in a reference data set taken 
with POLIS (right panel). Left panel, left 
to right: intensity spectrum, Stokes V of 
the destretched, speckle-deconvolved, and 
MOMFBD data, respectively. The display 
range of V is ±0.5%. Middle panel: the 
same for the second cut. Right panel, left 
to right: I and V (±0.1%) from the POLIS 
data, V plus 0.2% rms noise with a dis- 
play threshold of ±0.5%. The white mark- 
ers at the right border of Stokes V denote 
Vmax > 2.5<x, the red markers V max above 
the finally chosen significance level. 



These fluctuations of the polarity continue in the areas 
with weak or no polarization signals, but the three data re- 
duction methods show a slightly different behavior there. 
The destretched and MOMFBD spectra preferentially show a 
salt' n' pepper pattern with pixel-to-pixel variations typically for 
random noise (Fig. [9] and the magnification in Fig. [TO]). The 
speckle-deconvolved spectra on the other hand in some cases 
exhibit a very regular pattern of repeatedly alternating polarities 
extending over patches of 2 - 3 pixels each, for instance inside 
the white rectangle in the polarity of the speckle-deconvolved 
spectra in Fig. [TO] This regular pattern resem bles a few row s 
or columns from a "chequerboard" (Bello Gonzalez et al. 2009), 
but its regularity is rather surprising in the dynamic QS environ- 
ment. Because the pattern is less prominent in the MOMFBD 
data and almost completely absent in the destretched data, we 
tried to trace down its origin using the spectra along the spatial 
cuts marked by the vertical white dash-dotted lines in Fig. [9] 

The Stokes / and V spectra along these cuts are displayed 
in Fig. \TT\ Comparing the V profiles of the destretched data 
(second column in the left two panels) with those of both de- 
convolved data sets (third and fourth column), one can detect 
again a similar black-and-white chequerboard pattern in the de- 
convolved spectra, but this time the x-axis corresponds to the 
wavelength, not a spatial dimension. Most of these features that 
reach and exceed the display threshold of ±0.5 % amplitude are 
seen on 2 - 3 wavelength points, whereas the line width is larger 
than five wavelength steps because of the thermal and Doppler 
broadening, and over 2-3 pixels in the spatial dimension. The 
pattern therefore is generated by the reconstruction of the indi- 
vidual wavelength points and is not caused by any solar polar- 
ization signal. 

Its artificial origin is demonstrated also by the comparison 
to the right panel of Fig. \TT\ which shows a cut of identical 
length through the low-noise POLIS data set (see Sect. 12. 31) . 
The second column of the rightmost panel of Fig. \TT\ shows 
the original Stokes V spectrum with a threshold of ±0.1%. The 
Stokes V spectrum shows a polarization signal all along the slit, 
where all weak signals still have a larger spectral extent than the 
black/white pattern in the deconvolved Gottingen FPI spectra. 
The rightmost column of Fig. \TT\ displays the same spectrum af- 
ter an addition of a rms noise of 0.2%, which is typical for the 
Gottingen FPI data, with the same display range of ±0.5% as for 
the Gottingen FPI data. The added noise pushes about 40 % of 
the previously significant signals below the noise level, as for in- 
stance those inside the white rectangle near the middle of the cut. 
The clearly different appearance of the noise in the deconvolved 



Table 4. Final significance threshold in percent (top row) and 
relative to the rms noise cr in Stokes V (middle row). Area cov- 
erage of polarization signals above the threshold (bottom row). 



Speckle MOMFBD Destret- POLIS POLIS 
deconv. ching +noise 



% 


1.1 


cr x 


3.5 


of area 


18 



1.0 

5 
16 



0.8 
4.5 
13 



0.11 

3.5 
73 



0.6 

3 
35 (13 1 ) 



(1) after extrapolation to the spatial sampling of the Gottingen FPI 



spectra and the rightmost Stokes V spectrum suggests that the 
noise in the deconvolved spectra should not be treated as being 
Gaussian. It shows a spectral and spatial coherence over a few 
pixels in both dimensions caused by the noise filtering and the 
reconstruction process, respectively. The noise in the destretched 
data deviates from a Gaussian noise as well because of the spec- 
tral filtering, but rarely exceeds the display threshold of 0.5 % 
with noise peaks. 

We then determined suitable significance thresholds for the 
various data sets using the Stokes V spectra of Fig.[TT]as a guide. 
Profiles that exceeded the initially chosen 2.5 cr level are indi- 
cated by the white markers at the right border of all Stokes V 
spectra. This clearly selects many locations without any solar 
polarization signal in the spectra. The red markers correspond 
to locations that remain with the manually selected thresholds 
listed in Table [4] These signals can mostly be confirmed to be 
significant by eye as well. 

The rightmost columns of Figs. [9] and [TO] show the polar- 
ity of the locations of significant polarization signals inside the 
FOV of the Gottingen FPI when the rejection thresholds are ap- 
plied. They are now mainly limited to the signals that can be 
clearly identified in the map of the Stokes V amplitude (left- 
most column). The chequerboard pattern has disappeared almost 
completely from both deconvolved data sets. The fraction of the 
FOV that shows significant polarization signal is about 15 % for 
the Gottingen FPI data, in comparison to about 70 (30)% for the 
POLIS data (+noise) (bottom row of Table©. 

Stokes V amplitudes The reconstruction improves the spatial 
resolution and also amplifies both the polarization amplitudes 
and the noise level. Scatter plots of only the significant polar- 
ization signals showed an increase of the polarization amplitude 
relative to the destretched data by 20 % for the MOM FBD and 
40 % for the speckle deconvolution (cf. also Fig. IA.4I) . The in- 
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Fig. 12. NCP and intensity maps of the data sets. Bottom row, 
left to right, maps of the NCP for the destretched, MOMFBD, 
and speckle-deconvolved Gottingen FPI spectra, and continuum 
intensity. Top left: continuum intensity (top panel) and NCP for 
the POLIS data set. Top middle: same for the Hinode/SP data set. 
Top right: continuum intensity and NCP for the MHD simulation 
in full (left) and degraded resolution (right). All tick marks are 
in arcsec. 



crease of the polarization amplitude implies that the magnetic 
structures are better resolved in the deconvolved data sets with 
a higher filling factor inside a pixel. The amount of increase is 
in rough agreement with the ratio of the spatial resolution esti- 
mated in Sect. EJ (074/076 = 2/3, i.e., a 33 % smaller effective 
area per pixel for the deconvolved spectra). 

Net circular polarization The net circular polarization (NCP) is 
an important diagnostic quantity because it indicates gradients 
alon g the LOS in the magnetic field and the velocity field (see, 
e.g.. ISanchez Almeida & Liteslll992h . We defined the NCP as 
J VdA ignoring the polarity. Figure [12] compares the NCP maps 
of the various data sets: the three reductions of the Gottingen 
FPI spectra (bottom row), the POLIS and Hinode/SP observa- 
tions, and the MHD simulation in full resolution and spatially 
degraded to Of! 52 resolution (top row of Fig. [12]). The general 
shape of structures in the NCP of the slit- spectrograph obser- 
vations (POLIS, Hinode/SP) matches that in the MHD simula- 
tion with a spatial pattern on granular scales. For the destretched 
Gottingen FPI spectra, the NCP maps are rather noisy, but com- 
pare well to the Hinode/SP data set when noise of the level of 
the Gottingen FPI spectra would be added to the latter, showing 
connected patches of non-zero NCP values. 

The NCP map of the destretched spectra shows pixel-to-pixel 
variations outside the large-scale connected polarization patches, 
whereas both the MOMFBD and speckle-deconvolved spectra 
show a distinct background pattern of individual circles of a 
few pixel extent, more clearly seen for the speckle-deconvolved 
spectra. The noise level in the spectra affects the NCP strongly, 
such thus for instance in the left white rectangle of Fig.[T2lthe 
values from the three data reduction methods contradict each 
other: the destretched data suggest a negative NCP (black), in the 



- MHD Sim. 

-POLIS 

-SP 

■ Destretching 
■ MOMFBD 

■ Speckle deconv 




The observations are neither simultaneous nor co-spatial. 



Fig. 13. Histograms of the NCP values in the destretched 
data (thick blue), MOMFBD data (thick red), and speckle- 
deconvolved data (thick black). For comparison, histograms for 
Hinode/SP spectra (thin purple), POLIS spectra (thin orange), 
and for spectra from a MHD simulation (thin turquoise) are also 
shown. 



MOMFBD spectra the signal is basically noise, and the speckle- 
deconvolved spectra exhibit a positive NCP (white) there. The 
general noise level in the NCP of the speckle-deconvolved spec- 
tra is surprisingly smaller by an order of magnitude compared 
to the other two data reduction methods, even if the noise in 
the NCP should scale with the noise level in the spectra. The 
histograms of the NCP values in Fig. [13] show that the distri- 
bution of the NCP in the speckle-deconvolved spectra (thick 
black) matches better to all slit- spectrograph observations and 
the MHD simulation, whereas the destretched and MOMFBD 
spectra (blue and red curves) have much broader distributions. 
On closer inspection, the borders of the NCP patches in the 
speckle-deconvolved data seem to be artificially sharp (e.g., in 
the right white rectangle), changing from about zero to the dis- 
play threshold from one pixel to the next. The same behavior 
of the NCP also resulted for those MOMFBD test runs where 
the NB channel was not fully weighted in the deconvolution (see 
Fig.lA3l). 

The reason for this behavior of the NCP in some of the de- 
convolutions is not clear at once. The NCP is determined from 
the addition of several wavelength points, making it strongly 
sensitive to individual noise peaks. The shape of the distribu- 
tion (Fig. [13]), however, suggests that the noise in the NCP of, 
e.g., the speckle-deconvolved data is closer to that of observa- 
tions with higher S/N ratio than that of the other two data reduc- 
tion methods, even if the rms variation in the spectral dimension 
is higher for the speckle-deconvolved data. Adding up only the 
polarization signal in the blue or red wing of the line yielded 
nearly identical images for all data reduction methods, the sharp 
boundaries only appeared as soon as the two images were sub- 
tracted. We investigated the spatial Fourier power of the NCP 
maps, but this unfortunately also did not provide a clear clue to 
the observed behavior of the NCP. 



Cumulative polarization fraction Another method to compare 
data from different sources is the cumulative fraction of all pro- 
files that exceed a given polarization limit pn m . This fraction is 
a rather "robust" quantity because it is based on good statistics, 
and at the same time allows directly to see the detection li mit set 
by the noise level (see also Fig. 5 in Beck & Rezaei 2009). The 
noise level in the polarization degree can be read off from the 
location, where the curves start to deviate from 100%. Figure [T4l 
shows the result for the Gottingen FPI, POLIS, and Hinode/SP 



12 



K.G.Puschmann & C.Beck: Application of speckle and MOMFBD/MFBD techniques on 2D data 



50 % 150 % 
100 




i Thresh.; ivolution 



— ■ — ■ Simulation 

— — Simulation, degraded 

— 0.64 s GFPI SPECKLE 

0.64 s GFPI MOMFBD 

4.8 s, SP 

POLIS, 26 s 



Threshold MOMFBD p |im [% |J 

Fig. 14. Cumulative fraction of profiles with a polarization de- 
gree above p\ im . Thin black: POLIS data with 26 sees inte- 
gration time. Purple thin dash-dotted: Hinode/SR Blue dot- 
ted: MOMFBD Gottingen FPI data. Blue long dashed: speckle- 
deconvolved Gottingen FPI data. Thick red dash-dotted: MHD 
simulation. Thick red dashed: spatially degraded MHD simula- 
tion. The vertical red dotted line denotes the 50% fraction of the 
MHD simulation. The vertical blue dotted lines mark the sig- 
nificance threshold for the MOMFBD and speckle-deconvolved 
spectra. 

observations, as well as the original MHD simulation and the 
MHD simulation after degradation to a resolution of about Off 52. 

The respective curves and noise values generally follow the 
trend predicted by their integration time t, which is defined here 
as the total exposure time that entered into the measurement 
of the Stokes vector at one wavelength step (Gottingen FPI) 
or one slit position (POLIS, Hinode/SP). In comparison with 
the original and degraded MHD simulation (thick red curves 
in Figure IT4b . only the POLIS observation with t = 26 sees 
is able to measure the bulk of the polarization signals in the 
MHD simulation. For the Hinode/SP data (t = 4.8 sees), only 
the stronger half of the polarization signals of the MHD simu- 
lation (red vertical dotted line in Figure [14] at the 50 % level) 
are just at or slightly below the noise level. For the speckle- 
deconvolved/MOMFBD Gottingen FPI data (t = 0.64 sees), the 
cumulative fraction reduces to values below 100 % for a pu m of 
about 0.5 % of I c . Thus, only 30 % of the polarization signals 
present in the MHD simulations are above the noise level and 
can be detected. Only about 8 % (10 %) of the polarization sig- 
nals in the MHD simulation are above the significance thresh- 
old needed for the speckle-deconvolved (MOMFBD) Gottingen 
FPI spectra (vertical blue dotted lines in Figure [14]). The signals 
in the simulation, however, correspond to the comparably weak 
magnetic fields on granular scales. The Gottingen FPI and the 
Hinode/SP data show more locations with large polarization sig- 
nals (pu m > 1%) than the POLIS data, which results from their 
higher spatial resolution: for unresolved magnetic fields inside 
a pixel, the polarization degree increases with spatial resolution 
because it measures the fraction of polarized and unpolarized 
photons. 

5. Summary and Discussion 

We have applied two different image reconstruction techniques 
(speckle reconstruction/deconvolution, (MO)MFBD) to imaging 
and 2D spectropolarimetric data. For the imaging data (blue con- 
tinuum at 431.3 nm and Call H), both reconstruction methods 
yielded similar results with nearly identical rms contrasts, with a 
slightly better performance of the MFBD approach in the case of 



Call H imaging data with a low light level. For the reconstruc- 
tion of both the BB data (at 630 nm) or the NB data (Fe i line at 
630.25 nm) taken with the Gottingen FPI, the results differed sig- 
nificantly. The speckle approaches yielded a significantly higher 
continuum rms contrast (6.3 %) than the MOMFBD with equal 
weight for BB and NB channel (4.3 %), with a similar spatial 
resolution for both methods. 

A series of test deconvolutions (Appendix© showed a clear 
trade-off between contrast and noise level in Stokes V, where 
the solution of a MOMFBD with equal weights for BB and 
NB channel turned out to be the one with the lowest noise, 
whereas the speckle-deconvolution corresponds to one with 
medium noise and medium contrast. For a derivation of the mag- 
netic properties of the solar photosphere, the solution with low- 
est noise is to be preferred because the intensity contrast will 
enter stronger into the spatial variation of the derived tempera- 
tures than into the magnetic field properties. One would, how- 
ever, expect more similar results for the speckle deconvolution 
and the MOMFBD because the methods were applied to the 
same data set, although the final outcome depends strongly on 
the settings during the deconvolution (Appendix lAb. One rea- 
son for the difference between the speckle deconvolution and 
the MOMFBD intrinsic to the data used could be (small) dif- 
ferences in the optical path between BB and NB channel in the 
Gottingen FPI that the MOMFBD reconstruction is more sensi- 
tive to than the speckle technique (M. Lofdahl, personal note). 
Another reason could be uncompensated high-order ab errations 
dRouppe van der Voortet al. 2004: lScharmer et alJ2010h that in- 
fluence the contrast more than the spatial resolution. 

Comparing the deconvolved intensity spectra with a data 
set that was only corrected by applying a destretching algo- 
rithm, both deconvolution methods increased the spatial resolu- 
tion without introducing prominent data artifacts in the intensity 
profiles. 'Artifacts" in this context would be all patterns whose 
properties are not compatible with a solar origin, such as pixel- 
to-pixel variations above the limit given by the noise level (be- 
cause at the spatial resolution of the data solar features should 
cover consistently a few pixels), and small or large-scale patterns 
in the deconvolved spectra other than those seen also in the de- 
stretched spectra (e.g., spatially extended solar oscillations have 
to show up also without the deconvolution). Such cases are not 
seen in the line parameters shown in Fig. [8] The two parameters 
with a pronounced difference between destretched and decon- 
volved spectra are the line width and, to a lesser extent, the line 
asymmetry. The reason why these quantities change more than, 
for instance, the line-core intensity, presumably is that both are 
non-linear measures, where one additionally also expects a large 
intrinsic difference of the values inside granules or intergranular 
lanes. The increase of the spatial resolution by the deconvolu- 
tion therefore has a profound effect on these two quantities, but 
the change seems to be of solar origin. The grainy pattern in, 
e.g., the line width of the speckle-deconvolved data is so weak 
that it does not qualify as significant because its amplitude will 
presumably not be noticeable in an analysis of the data. 

An investigation of the polarization signal revealed, however, 
one of the drawbacks of the deconvolution, i.e., the amplification 
of the noise level. All of the spectra have to be noise filtered in 
the spectral domain to be useful. For the speckle-deconvolved 
data, spurious polarization signals of up to 0.5-1.0% of I c re- 
mained even after a rather strong filtering. The high noise level 
in the deconvolved data can affect the analysis of the polarization 
signal producing for example a chequerboard pattern in the po- 
larity of the magnetic fields that is absent in the destretched data. 
The pattern could be traced back to individual noise peaks in the 
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Stokes V spectra that are not caused by a solar polarization sig- 
nal. From a visual inspection of the spectra and a comparison to 
a reference slit- spectrograph spectrum with a noise level lower 
by one order of magnitude, we determined significance thresh- 
olds for the deconvolved data. 

For the speckle deconvolution (MOMFBD), signals with 
a polarization amplitude below about 1.1% (1.0%) of I c can- 
not be taken as significant. Application of the correspond- 
ing thresholds removed the chequerboard pattern from the 
data. The 1 % -level needed for the deconvolved Gottingen 
FPI spectra prevents the detection of weak and diffuse mag- 
netic fields. Using the terminology of Stenflo (2010) that so- 
lar magnetic fields in the quiet Sun exist in either a "col- 
lapsed" (partly evacuated and concentrated magnetic flux) or 
an "uncollapsed" (diffuse and weak magnetic fields) state, the 
Gottingen FPI can only sense the collapsed fields. This might 
explain why previous studies of the QS with the Gottingen 
FPI (e.g.. iDommguez Cerdena et alJl2003l 12006) found a much 
larger fraction of magnetic fields in the kG range in the 
QS than slit- spectrograph observations of visible or infrared 
spectral lines (e.g., [K homenko et al.l [200 3. 2005; Rezaei et al 



2007; Orozco Suarez et al. 2007; Martinez Gonzalez et al. 200L 
Beck& Rezaei 120091) . Using a MHD simulation with weak in 
ternetwork magnetic fields as a reference, a noise level below 
about 0.25% of I c is needed to detect the weaker half of the po- 
larization signals (for the case of IBIS cf. Wo ger et all 2009). 
This level is not met by neither the Gottingen FPI data nor an 
observation with the Hinode/SP with 4.8 sees integration time, 
it requires significantly longer integration times to achieve a suf- 
ficiently high S/N ratio. 

The distribution of NCP values for the speckle-deconvolved 
data matches that of the observations with higher S/N better than 
the other two data reduction methods but also has seemingly ar- 
tificially sharp boundaries, a feature that showed up in the NCP 
maps of all deconvolutions where the information from the NB 
channel was not or was only partially used. The behavior of the 
NCP should therefore be related somehow to a step that is spe- 
cific to these deconvolutions, but neither the MOMFBD weight- 
ing both BB and NB channel equally nor the destretching of the 
spectra. One candidate would be the way of determining and ap- 
plying the OTF. In the destretching of the spectra, the OTF is not 
applied directly to the spectra, they are only spatially aligned to 
the reference BB image. In the MOMFBD weighting both the 
NB and BB channel equally, the OTF is derived using the full 
information from the NB channel, whereas in all other cases the 
information from the BB channel is dominating in the deriva- 
tion of the OTF. In the presence of (minor) optical differences in 
the light path of the BB and NB channel, the OTF based on the 
BB channel might therefore not fit perfectly to the NB channel. 
The application of the OTF to deconvolve the NB spectra could 
then result in a small additive effect that only shows up when 
wavelength-integrated quantities such as the NCP are calculated 
from the spectra. The NCP maps of deconvolved spectra with a 
reduced noise level in the NCP - contrary to their increased rms 
noise in the Stokes V signal relative to the destretched spectra - 
and the overly sharp boundaries in the spatial domain fulfill the 
definition of showing "artifacts" as given above. 

6. Conclusions 

To determine the "best" approach for the deconvolution of po- 
larimetric spectra, one has to choose a compromise between spa- 
tial resolution, intensity contrast, and the noise level in the po- 
larization signal. The destretched spectra provide a reference for 



the lowest achievable noise level that in our case is roughly pre- 
served by the MOMFBD of the spectra with equal weight for BB 
and NB channel. For the future usage of the Gottingen FPI for 
imaging spectropolarimetry it, however, seems to be necessary 
to improve the S/N ratio, to be able to lower the final significance 
thresholds. Especially for weak spectral lines with a small line 
depth, and hence, small polarization amplitudes, or strong chro- 
mospheric lines with a low line-core intensity, the most common 
polarization amplitudes are far below the 1 % level and thus can- 
not be detected at a low S/N. 
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Appendix A: Deconvolution test runs 

We run a series of deconvolutions with different settings. For the 
speckle deconvolution, we used on the one hand no Fourier filter- 
ing during the deconvolution and on the other hand an optimum 
filter (no. 1 and 2 in Table lA.li) . A series of MOMFBDs was done 
starting from one weighting only the BB channel in the deconvo- 
lution (no. 3). The spectra of no. 3 were then Fourier- filtered post 
facto with a Hamming filter in the Fourier power of individual 
wavelengths and Stokes parameters (no. 4 to 9). The post- facto 
filter was set to cu t all power beyond a given spatial scale (second 
row of Table lAJl) . Instead of a post- facto filter, we also changed 
the noise threshold in the MOMFBD already (no. 10 and 11), 
still only weighting the BB channel. All other MOMFBDs used 
the default noise filter (NF=1). In the final series, we included 
the information from both the BB and NB channel with an in- 
creasing weight for the NB channel from 0.1 to 1 (no. 12 to 15). 
The last entry (no. 16 in Table lA.li) corresponds to the destretch- 
ing of the spectra. The spectral noise filter of Sect. 13.41 was not 
applied, the values of the noise in V/I at continuum wavelengths 
(630.275 nm to 630.291 nm) thus correspond to the first column 
of Tabled] 
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Fig. A.l. Continuum intensity {bottom row) and average unsigned Stokes V/I signal < 1.5% (top row) for the deconvolution test 
runs. 

Table A.l. Intensity contrast of the NB channel and rms noise in V/I in the test runs. 
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method 
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weight NB 
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filter 
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NF2 
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contrast [%] 


11.1 


8.5 


5.8 


4.9 


4.9 


4.3 


3.6 


rms V/I [%] 


0.89 


0.59 


0.54 


0.41 


0.41 


0.32 


0.29 



A.1. Reasonable parameter range 

The destretched spectra (no. 16) provide the minimum noise 
level in V/I (0.29%) and the minimum contrast (3.6%). The 
maximum contrast is provided by the speckle deconvolution 
without filter (no. 1, 6.3 %) and the MOMFBD deconvolutions 
weighting only the BB channel without additional filter or with 
a slightly enhanced noise filter (no. 3 and 10, 10-11 %). The 
"best" settings have to yield contrast values within this range, 
without a rms noise in V/I below that of the destretched spec- 
tra. The upper limit of the noise level is not restricted (> 1 % for 
no. 3). 



A.2. Contrast vs. noise 



Figure lAJl shows a subsection of the full FOV in continuum in- 
tensity (bottom row) and unsigned average Stokes V/I signal (top 
row) for the different deconvolution settings. The variation of the 
noise level can be clearly seen. Variations of the intensity con- 
trast are hidden because each panel is scaled individually, only 
the cases with a degraded spatial resolution (no. 8,9, and 16) can 
be distinguished at once. Figure \AJ2\ visualizes the trade-off be- 
tween noise and contrast for the case of the post- facto filtering in 
the MOMFBD spectra. The noise level falls much more rapidly 
with the increasing spatial scale at which the power is cut than 
the intensity contrast. The noise level in the destretched spectra 
sets an upper limit of about 07 7 for the filter cutoff. The con- 
tinuum intensity images in Fig. IA. II set a stricter limit. Between 
no. 7 and no. 8 (0735 and a/45), a degradation of the image qual- 
ity can be clearly seen, whereas for the previous steps (no. 3 to 
7) the degradation is not prominent. Larger filter widths signif- 
icantly smear out the spatial resolution or fall below the noise 
level of the destretch (no. 9). 

All MOMFBDs where the NB channel was not fully 
weighted give a slightly more diffuse image, which is more 
prominent in images of the full FOV (not shown here). All 
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Fig.A.2. Dependence of intensity contrast (bottom) and rms 
noise in V/I (top) in the MOMFBD spectra on the cutoff scale of 
a post- facto Fourier filtering. The horizontal dotted lines give the 
corresponding values for the speckle-deconvolved, MOMFBD 
(no. 15), and destretched spectra, respectively. 



MOMFBDs with only a partial weight for the NB channel 
(no. 12 to 14) had some more or less pronounced artifacts in the 
spectra, mainly some individual badly reconstructed subfields. 
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Speckle deconv. 
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Fig. A.3. NCP maps of some of the test deconvolutions. The dis- 
play range is ± 0.15 pm. The FOV corresponds to the one used 
in Fig. El 
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Fig. A.4. Stokes V/I(A = 630.245 nm) amplitude of deconvolu- 
tion no. 7 vs. that of no. 15. The short-dashed line denotes a 
unity slope, the long-dashed line a slope of 1.5. 

A.3. Net circular polarization (NCP) 

The NCP is a sensitive measure to determine the spectral pu- 
rity of the deconvolved spectra, in addition to the noise level. 
The scatter in the NCP values should increase with the increas- 
ing noise level in V/I. It, however, turned out that this is not 
necessarily the case. Figure lA3l shows a subsection of the NCP 
maps of the full FOV for some deconvolutions. For the methods 
no. 14 (MOMFBD weighting NB with 0.5), no. 11 (MOMFBD 
only weighting BB, NF = 2), no. 7 (MOMFBD only weighting 
BB, post- facto filter cutting at 0745), and no. 2 (speckle de- 
convolution) the scatter in the NCP is strongly suppressed and 
does not scale with the noise level. The NCP shows artifically 
sharp boundaries in the NCP maps, changing from about zero 
to the display threshold of ±0.15 pm from one pixel to the next, 
e.g., inside the white rectangle marked in Fig. IA.3l This behav- 
ior is especially seen for the deconvolutions where the infor- 
mation of the NB channel was not or only partially used, with 
the possible exception of no. 7. Neither the destretching nor the 
MOMFBD weighting both NB and BB channel equally (no. 15) 
show the sharp boundaries. The latter deconvolution is the only 
one, where the noise in the NCP of the deconvolved spectra in- 
creases with the noise level in V/I (compare the first two maps 
of Fig. DOl) . The spectral noise filter of Sect. 13. 41 was applied be- 
fore the calculation of the NCP, but it has no effect on the NCP, 
and especially does not cause the artifical behavior of the NCP. 
Data without the spectral filter show the same patterns. 



the largest fraction of pixels inside the FOV, the relation gives a 
unity slope between the two co-spatial polarization amplitudes. 
For all high-amplitude signals (\V/I(A = 630.245 nm)| > 5 %) 
the slope changes to 1.5, implying a relative increase of the po- 
larization by 50%, e.g., 7.5% instead of 5 %. A direct compar- 
ison of the two V/I (A = 630.245 nm) maps showed that the po- 
larization amplitude in the central part of the network patches 
is enhanced in no. 7 and reduced in their surroundings, i.e., the 
"halo" produced by spatial smearing was better removed than in 
no. 15. The increase of the polarization amplitudes is, however, 
the only argument in favor of the deconvolutions with a contrast 
above that of no. 15. 

A.5. Final choice of settings for MOMFBD 

The results of the various test runs differ mainly in two parame- 
ters, the rms contrast in intensity and the noise level in V/I, with 
a clear trade-off between the two parameters. The spatial resolu- 
tion seems to be virtually unaffected as long as no strong spatial 
filtering is involved. Because finally the polarimetric properties 
of the spectra are the observed quantities that are used to derive 
the properties of the solar atmosphere, we chose the MOMFBD 
with the lowest noise level (no. 15) in which NB and BB channel 
entered with equal weight in the MOMFBD. Besides the noise 
level, this choice also yields more reasonable values of the NCP 
than some of the other deconvolutions. The impact of the choice 
on derived solar properties can be easily tested by, e.g., compar- 
ing the results of an inversion of the speckle deconvolved spec- 
tra (medium noise and contrast) with those of an inversion of the 
spectra of MOMFBDs no. 3 and 15, and the destretched spectra. 
The rms contrast should enter into the spatial variation of the 
temperature with hopefully a minor change of the derived mag- 
netic field properties. This is, however, beyond the scope of the 
present investigation. 



A.4. Amplification of the polarization amplitude 

The deconvolution does not only scale up the noise level with 
increasing intensity cont rast, but also the amplitude of the po- 
larization signals. Figure \AA\ shows a scatter plot of the polar- 
ization amplitudes at 630.245 nm in method no. 15 (NB fully 
weighted) and no. 7 (only BB weighted, post- facto filter cutting 
at 0745). The method with the larger contrast and noise (no. 7) 
has increased polarization amplitudes, but not on all pixels. For 
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